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1. INTRODUCTION AND SUMMARY 

A special materials processing furnace is being developed by NASA- 
MSFC for use in forthcoming Spacelab missions to study the solidification 
under closely controlled conditions of various sample materials in the absence 
of gravity. The samples are to be rod shaped and will be subjected to both 
heating and cooling simultaneously with the heating and cooling distributed 
along the length of the rod by annular heat pipes concentric with the sample 
rod. The heat pipes will serve to maintain near-uniform temperature dis- 
tributions along the hot and cold sections of the rod with a very steep tem- 
perature gradient in between. 

The thermal model is based on a Lockheed-Huntsville developed Thermal 
Analyzer computer program. The model was developed to be very general to 
enable the simulation of variations in the furnace design and, hence, serve as 
an aid in finalizing the design. 

The thermal model is described in Section 2, with a user's guide given 
in Section 3. Some preliminary results obtained in testing the model are 
given in Section 4. 
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Z. THERMAL MODEL 

The thermal model for the furnace was based on a number of preliminary 
design drawings and sketches provided by NASA-MSFC. A simplified schematic 
of the furnace is shown in Fig. 1. We concentrated our detailed thermal model- 
ing efforts in the inner core region, including the sample, heat pipes, core tube 
(heating elements) and heat sink. Some of the outer region was not modeled 
exactly, since this would introduce considerable complexity into the model 
without significantly enhancing accuracy. 

In modeling the furnace, we divided the overall model into three regions. 
The center section contains the core tube, hot heat pipe and surrounding insu- 
lation and housing and is labeled as Region 1. The end section opposite to the 
cold heat pipe and heat sink contains primarily structural stiffeners and insu- 
lation and is labeled Region Z. The end sectio l containing the cold heat pipe 
and heat sink is labeled Region 3. 

Heat is generated in the core tube and is distributed to the sample rod 
after being transferred through the muffle and hot heat pipe. Some heat is 
lost by conduction through the multilayer insulation to the external housing, 
where it is collected in a cooling jacket. The heat flowing through the sample 
into the cold heat pipe and heat sink is collected and removed through the heat 
sink. 


A detailed schematic of the Region 1 nodal network is given in Fig. Z 
with a table of conductors. In this coniiguration, the sample rod is fully 
extended into the heated region. The position of the sample can be altered 
to simulate translation from the heated region into the cc led region. We 
assumed that the heat pipe is an infinite conductor and, hence, is isothermal. 
The region within the muffle is partially pressurized by helium to enhance 
heat transfer by conduction. We modeled heat transfer in the radial direction 
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Fig. 2 - Schematic of Region 1 Nodal Network and Table of Conductors 
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within the muffle by a combination of conduction through the helium gas 
medium and radiation acrois the gaps separating the sample, heat pipe and 
muffle. Outside ox the muffle, we modeled heat transfer in the radial di- 
rection by radiation only. Radiation through the multilayer foil is given by 

■> ■ TT^TTn (T t - T 2> "» 

wheic q is the radiant heat flux, e is the emissivity of the foil, n is the num- 
ber of foil layers and and T, are the temperatures on the surfaces bound- 
ing the n foil layers. 

Heat transfer in the axial direction was modeled by conduction along 
the structural elements. 

A detailed schematic of the Region 2 nodal network and conductors is 
given in Fig. 3. The configuration of stiffeners and multilayer insulation 
which we modeled in Region 2 is a simplification of the actual design shown 
in Fig. 1. This simplification was made to facilitate the modeling of heat 
transfer in this region without unduly sacrificing accuracy. Heat transfer 
was modeled by conduction along the stiffeners and radiation through the 
multilayer insulation. 

The Region 3 nodal network and conductors is shown in Fig. 4. The 
sample is thermally connected to the cold heat pipe only if it is translated 
into the cooled region. As in Region 1, heat transfer within the muffle is 
by a combination of conduction through the helium gas and radiation be- 
tween structural members. Outside the muffle, radial heat transfer is by 
radiation only, and axial heat transfer is by conduction along structural mem- 
bers. The cold heat pipe, like the hot heat pipe, is assumed isothermal. 

In Regions 1, 2 and 3, the external housing is maintained at a constant 
temperature by a cooling jacket. The furnace is assumed to be controlled by 
maintaining the four core tube sections at fixed temperatures. The cold heat 
pipe is maintained at a fixed temperature by cooling from some arbitrary 
heat sink. 
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Figure 4 - Schematic of Region 3 Nodal Network and Table of Conductors 
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3. COMPUTER PROGRAM 

A complete lilting of the computer program for the thermal model is 
given in Appendix A. The program is coded in FORTRAN V language for 
the NASA-MSFC Univac 1108 computer The program can be run for various 
core tube temperatures and sample translations with various heat sink and 
cooling jacket temperatures by varying the program inputs. The thermal 
model can be changed to simulate furnace design changes by modifying the 
program. 

The main program contains the thermal analyzer coding, but calls sub- 
routines for defining initial conditions, nodal networks, capacitances and con- 
ductances. Initial conditions for the temperatures are computed in subroutine 
CONI. The initial temperatures in Region 1 are computed based on fixed core 
tube and housing temperatures assuming radial conduction at equilibrium. 

The initial temperatures interior to the core tube in Region 1 are made equal 
to the core tube temperatures. The temperatures in Regions 2 and 3 are 
initialized at the cooling jacket, heat sink cr some other temperature near a 
Regio:i 1 temperature. The idea is to start the program calculations at initial 
temperatures as near as possible to the actual equilibrium temperatures. This 
approach to initializing temperatures reduces the required program run time. 

The linking between nodes in the nodal network is defined in subroutine 
LINKUP, and conductances for each link are defined in subroutine COND. In 
COND, there are basically three kinds of conductances considered: (1) radia- 
tion, (2) solid conduction, and (3) helium conduction. The radiant conductance, 
CD, for a single gap is given by 

CD ■ - ■ 0(T ^ ^ (T ! + T 2* A 
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where f and , are the emiasivities and T and T , are the temperatures of 

1 b l C 

the two surfaces, o is the Stcphan-Boltzmann conatant, and A ia the area. 
The conductance for multilayer lnaulation ia given by the aame expreaaion 
divided by the number of foil layera. 

The conductance of aolid conductors ia given by 



where K ia conductivity, AX ia thickness and A is c roas -sectional area. The 
K values are read into the program in table form in subroutine DATA as a 
function of temperatui e. The tabular data is interpolated in subroutine INTP. 
Provision is made for 16 tables in the program giving conductivity and specific 
heat for stainless steel, alumina, tungsten, molybdenum and tantalum alloy. In 
addition, tables are set aside for the sample material and the hot and cold heat 
pipes. The tables are described in the listing of COND In Appendix A. 

Conductan >re also calculated for helium conduction in the region 
contained with; the muffle. These calculations are noted by comments in 
the rrogram listing. 

Since the program computes an equilibrium solution, the values of the 
capacitances are not significant to the final solution except when used to flag 
nodes where temperatures are held constant. When the value of the capaci- 
tance is -1, the node is flagged as a constant temperature node. This is done 
in subroutine CAP. Subroutine HTRT is provided for defining heat generation 
rates. In our case, however, we set all heat rates Q equal to zero and, in- 
stead, use constant temperature nodes at points of heat generation or loss. 

The time step in the numerical solution is recomputed to ensure sta- 
bility at each step in the solution. This is done in subroutine TSTP. 

The program inputs and their formats are given in Table 1. 
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The program output lists the program input data, the conductors and 
their nodal end points, the nodes and the connecting nodes and conductors, 
and the tables of material properties prior to entering the thermal analyser 
loop. Alter entering the thermal analyser loop, a listing is made of tem- 
peratures, heat rates and conductance values at specified intervals. After 
completion of a run, totals are made of heat transferred to the cooling jacket 
(QCASE) and to the cold heat pipe heat sink (QCOOL). The sum of these heat 
flow rates is the furnace requirement (POWER) The heat flow rates QCASE 
and QCOOl, an printed out in units of kcal/sec, and the power level POWER 
is printed out in kilowatts. 
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1 . 

4. SAMPLE PROBLEM 

A sample case was run for a tungsten sample rod translated 5 in. into 
the cooled region. The four core tube temperatures ere assumed to be 
8 38 C (2000 R), with the external housing and cold heat pipe cooled to 27 C 
(540 R). The input cards for this case are listed in Appendix B. The values 

f 

of PRIN and NEXCIT are not significant to our model and were arbitrarily set 
equal to 100. NPRl was set equal to 50 to print date every 50 iterations. 
TSCON was set equal to 2 to assure a stable solution. TMAX was set equal 
to 100, since we found by experiance tbit equilibrium is reached by that time. 

The calculated temperatures at equilibrium are plotted in Fig. 5 along 
a radius at the midpoint of the core tube. The calculated sample temperatures 
are plotted in Fig. 6 along the sample length. The total heat lost to the external 
hording cooling jacket QCASE was found to be 146 kccl/hr, and the heat flow 
to the cold heat pipe was found to be 255 kcal/hr. This corresponds to a total 
power level for the core tube of .47 kW. 
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Core Tube Temperatures = 838 C (200 R) 
Cold Heat Pipe Temperature = 27 C (540 R) 



Fig. 6 - Temperature Distr ibrtion along Tungsten Sample Rod 
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5. CONCLUDING REMARKS 

The thermal model developed during this effort accurately simulates 
heat transfer in the baseline furnace design, and is sufficiently general to 
model a broad range of design modifications. The model can be utilized to 
perform heat transfer analyses for various anticipated test conditions to 
determine required furnace power levels, cooling rates and sample tem- 
perature distributions. 

It is recommended that experimental data be obtained on the conductance 
of the multilayer foil insulation. This updated data can be incorporated into 
the present thermal model with little difficulty. 
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Appendix A 

LISTING OF COMPUTER PROGRAM FOR SPECIAL 
MATERIALS PROCESSING FURNACE THERMAL MODEL 


LOCKHEEO • HUNTSVILLE RESEARCH & ENGINEERING CENTER 



n r> 


LMSC-HREC TR D697553 


C* • • ***** 

COMMON/THANS/ li T K ANS 

CuMMON/XiNlT/ ICONl »CDb7.CU7b.OUto*.CUV10 

T H fc W M A L. ANALYSE* tJY JUHN MjNU 


U l MENb I UN T ( 2 JO ) »C ( 2 00) .CO ( 500 > »U<200 > * TO < 200 ) 

DIMENSION L. ( bOO • 2 ) »M( 200* I 5) 

DIMENSION TC0i<E(4) 

70 1 FORMAT (215. El 0.0 ) 

7j£ FORMAT (415) 

703 FJRMAT(E1''.0»2I b»4E10.0) 

704 FORMAT ( 4E 1 0.0 ) 

70 5 FORMAT ( 1 HI . t PU 1M« • *F 1 0.0 «2X. « NPR 1 > ' * 15«2X. • Nfc>CT« • • 1 5 • £X 

* FlC. 1 .2X. *TMAX*-*F 10.0) 

70 6 FORMAT ( 1HO. t I CAsE* • «F I0.0»2X « • ICCOL* • *F 10.0»2X.4 ( • TC0«E ( 
■1* F 10.0.2X. ) ) 

Vu 7 FORMAT ( IMO. • N 1 • . 12) 

1 lb.JX.bHNl T * .15) 

* F | 0 • 0 » 2X « ) ) 

To 7 FORMAT ( 1 HO * * M f KAN5» • « 12) 


1 





15. 

3X.5HNIT s 

71 1 FORMATE 1H 

* ( 2HT ( 

9 

1 5»2H) = 

» 2X 

•E10.4.3X. 

1 


2HT ( 

9 

15* 2H) a 

• 2X 

. E 10.4 «3X. 

2 


2HT ( 

» 

lb.2H)r 

. 2X 

. E 10.4. OX. 

3 


2HT ( 

9 

15. 2H)« 

• 2X 

.E10.4.3X. 

4 


2MT ( 

• 

15. 2H) = 

• 2X 

.E10.4.3X > 

712 f^rmat< 

1H 

* <2MC< 

• 

1 b. 2H) * 

*2X 

.FB.3.3X. 

1 


2HC( 

t 

15. 2H> = 

• 2X 

, Ftl. 3 » 3X . 

2 


2HC( 

t 

1 b.2H) = 

• EX 

.FB.3.3X. 

3 


2MC( 

9 

15* 2H) = 

»2X 

» Fb. 3 * 3X . 

4 


2HC( 

9 

1 5.2H) s 

• 2X 

»FB. 3*3X) ) 

TlJ F ORMAT ( 

1H 

« ( 211*. l 

• 

l 5 » 2 H ) = 

» 2X 

. L 1U.4.3X. 

1 


2HM 

• 

15. 2H) = 

. 2X 

. E 1 0 . 4 , JX , 

2 


2HK.I 

9 

!5«2H) = 

.^X 

•E10.4.3X. 

3 


2 HK. ( 

9 

1 5 . 2h ) s 

. 2X 

.E 10.4. JX, 

4 


2H<1 

9 

1 b.2H) c 

. 2X 

.E10.4.3X) 

7 14 F ORivAT ( 

1H 

* < 2HJ( 

9 

1 5.2H) * 

»2X 

. E10.4.3X. 

1 


2HJ« 

9 

I5.2H) * 

. 2X 

. E10.4.3X. 

2 


2HJ( 

9 

1 5 . 2 H ) = 

»2X 

» E 1 0. 4 » JX . 

3 


2H01 

9 

15. 2H)= 

. 2X 

. E 1 O . 4 , JA , 

4 


2H0( 

9 

15«2H> = 

* 2X 

. E 10.4. JX) 


15) 


7 1 5 F ORM A T ( 
71b F" ORM AT < 

71 7 FORMAT ( 

718 FORMAT ( 
1 

2 

3 

4 


1 HO « 3 1 b ) 

1 HO » 2 1 7 / ( 1H ,2015) ) 

1 Hi , 24HC0NDUCTOR NOOE- l NOOE-2) 
lHl,pbHNOL)E NUMBER OF RESISTORS.// 

1H , lOHNOOE — CON. * 5X , lOHNOOE — CON. »5X. 
1H «i 0HNU0E--CON. , 5X , lOHNOOE — CON. ,5X. 
1H ,j QHNOOE — CON. , 5X , lOHNOOE — CON. «5X. 
IH . i 0HN00E--CON. »bX. 1 OHNODE--CON. »bX> 


T bCON = • * 
1 1 » • ) * • * 
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720 FORMAT!//.', X, ldrlMA* TEMP CHANGE ■ » F l 0 • S • J* • 9HA T NOOE 

1 26Hi«OM OF Abb TtMP CHANGtS * »E15.6> 

7o0 FORMAT ( l rt • ( t>MJ0 JT ( . I ts » 2H > * » 2X • E 1 2* A • 3* • 

\ bMJOOT ( • 1S.2H) « .2X.E 12*4 .JX. 

2 *>HJ00H • lbti!H>»»*X*£lii*A»JX« 

j b MO 001 ( • | *_>» 2H )«»2 X»l12.4»3X) ) 

2 000 oUnT I Nut 

PtAOIS. 70 J ) P K I <M * NPR 1 «NEXC T .TSCON. r MAX » T L ALE » T COOL 

PEADl'S.704 ) ( T COKE 1 I > * I « 1 .4 ) 

PEAO! S. 702 InTKAiHS 

WW| TE I 6« 705 >PK IN* NPR J .NEXCT « TSCON. TMAX 

I T E 16 * 706 ) I v.Ao£ • Tcool • < 1 • TCOpe ( I ) . I ■ 1 . 4 ) 

WKJ TEI6. 707 jNTRANS 


IF ! PRIN.EU.O *0 ) OO TO 4000 
NNDS* l 60 
NCD6»347 
KKFt-G« 100 

CALL CONI ! T , TLOKt * TCASE* TCOOL .CO. L ) 
CALL LINKUP! L I 

C SORT OUT COnDUC T OPS AND AOUACENT NODES 

C 

c* **»##** 

OO 606 IN = I.NnOS 
NK ■ 0 

00 607 IP = 1 .NCOS 

IF! LI IR» 1 ) »nE • I H ) GO TO 60b 
NP = NK ♦ 1 

11 « NK *2 

12 » NP*2* 1 * 1 

Ml IN. II) = 1 ( IP ,2) 

Ml IN. 1 2> = IK 

GO TO 609 

605 CONTINUE 

IF|H IP. 2) .Nt. IN) GO TO 609 
NP s NK ♦ 1 

1 1 = NK*2 

12 * NP*2-M 
M1IN.II) = l(IK.I) 

M 1 I N . I 2 ) » IK 

6u9 CONTINUE 
6 u 7 CONTINUE 

M I J N . 1 ) at NR 

606 COnT I NUE 
C*«*****»**« 

c 

c pp i nt uut InIIIAl infupmation 
c 
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c • • • *•• • *** • 

C Tt ( C« 71 7 ) 

i>o ‘j 0«: I * i ♦ I'iCDb 

! If: (6* 7 I‘j) I »LC 1*1 > «L( I ••! > 

‘>0*- CONTI NUt 

WW| Tt I 6« 71 «> 

00 bOO 1 3 1 «nNl>s 
|N>.4< I • 1 ) *2<fl 

*W| T£ ( 6« 71 6 ) I * « M < I « J) » J«1 » IN) 
bOU COnTI NU t 
TIM « 0«0 
NlT ■ 0 
IPT* 1 
NPR * 0 
OT = 0*0 
NCTWL » 0 

MAxTNQ « 0 
TMAXCG * 0.0 
TSMCG » 0«0 

HEAD IN TAWJLAK data 

call data 


b f AWT LOOP 
OO 1012 1*1 ,N*Dh 
C( I ) =2.0 
1012 CONTINUE. 

1O00 CONTINUE 

CALL CAP(C.T) 

CALL COND< CD * L » T ) 

c 

c****»**» 

c 

C COMPUTE HLAT KATtb 

call HTWT < T , r M»OT «NNOb«U) 
lb2 CUNT INUt 
C 

c 

C PW|NT OUT HESOLTo 


1 E ( NPP 1 * ( NPR / i>|p K 1 ) • NE • NPR ) GO TO 160 
WR | TE (6*710) T I M »DT • NC TRL • N I T 
WPI TE l fc. 720 ) 1 MAXCG.MAXTNO* TbMCG 

»h| Tt(fc*71 1 ) ( i«» T IN) <N>1 .NNOb) 

*R 1 T E 1 6 • 7 1 4 ) liViOlN) . N = 1 • NNDb ) 

*R I TE < C * 7 I J) 1 it * CO * N ) tNs 1 * NCDb ) 
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too CUNT I NUt 

NPU * NPM ♦ 1 

too CUNT INUL 

c 

c* * 

c 

C COMPUTE time i» T tP 

CALC TblP<C.CU»l*CI*L»M.NNDb.NPW.PPlN»T|M. IPT • 1 bCON « ITT ) 


MAxTNU * 0*0 
Tiv.AXCG « o.o 
TSMCG * 0«0 
TIM a TIM ♦ o T 
NI T*N I T-M 
DO U00 1 a i ,\,\U5 

tcg » Au^t t< i )- r J< t ) ) 

1F( TCG.tT. TmAXCO) GO TO 94b 
MAyTND = I 
TMAXCG a TCG 
V4b CUNT l NUt 

T bMCG* TbMCG* T CG 
TU( I) » T ( I ) 

eoo cunt i N ut 


COMPUTE NfcW TtMPfcWATUPES 


OU 900 1 a ] «imNDS 
IF ( C ( 1 ) )90J, 902.901 

transient CALCULATION 

90 1 CUNT l NUt 

NADU a Ml J « i ) 

F 1 = 0 
Ft = 0 

DO 803 U a i ,„AUJ 
Nk | s J*2 
NP2 s U*2+l 
UN s M ( I « NW J ) 

JC = Ml I • NP? > 

Ft a Ft ♦ TOI UN) *CD'. JC) 

F2 a F2 ♦ Toll ) *CD< JC ) 
bo J CUNT I NUt 

T ( | ) x TO! I >-MF l-F2+3« I > >*DT/CI I ) 
Gu TO 900 


A -4 


LOCKHf EO HUNTSVILLt REST ARCH K ENGINEERING CENTER 


non 


LMSC-HREC TR D697553 


c 

C STfAUY STATp CALCULATION 

C 

90 2 CONT I NUL 

NAO J ■ M( | « | ) 

F I * 0 

Ftl * 0 

DO bOA J » 1 • •'■AD J 

NW1 ■ 0*2 

NK2 ■ J*2*l 

jN ■ Mil* N k I I 
JC « Mill NH2 ) 

FlaF | ♦TOt JN ) *CD( JC ) 

F2 a F2 ♦ CDtJC) 
bOA CONTINUE 

T( 1 l«(FI+J( | ) > / F 2 
9uJ CUNT INUt 
900 CONTlNUt 

MOT HE. A I PIPE. TtMPLWATuPL AS SUM I >'*Q iNFINITt CuNOUCToW 

jUM I a 0 • 

bOMc 1 a 0 • 
bUMJ«0t 
bUM 4 iO« 

1 MAX a 12-N*“MANb 
00 290 1 — 1*1 MAX 
K= t ♦ I 0 * < l - l ) 

SUM 1 a SUM 1 ♦ T ( K ) * CO » K > 

SUM2«SUM2+CD( N ) 

290 CONTINUE 

OO 291 1 3 I • 1 2 

K = 4 + 10*11-1) 

SUMO * SUM J ♦ T ( Nt 1 ) *COl < ) 

SUM A a SUM4 tCOU) 

29 « LONT INUC 

SUMS* SUM J + I ( 1 A 1 ) *CO< 3*7 ) 

SUM4*SUM4+COl 3 A 7 ) 

Tills! SUM 1 *SUM3 ) / ( SUM 2 + SUM 4 ) 

OO 300 I *2 • I 1 2 . 10 
T < | > a T ( 3 ) 

T ( | ♦ 1 l*T( Jl 
r ( r +2> »T < 3 > 

J00 CONTINUE 

JF< T IM.Ur. . TmAa ) GO TO 1000 
NPTSsNPhi 
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HtftI MAItb 

ucasc«o* 

DO 320 l«9. | lV.iO 
UCASE«UCAst+Co( l)«IT( ||-TCAi>t) 

JtC lOn! 1 Not 

UCA^b^CAI.t4(.bM 24»7)*( T( | iV)-1 < 120)1 
OCAhE-OLAhE-t-CDl 252>*< T( 130)-T( 121 ) 1 
OCAbE«ULASfc*Co< 2i>A >*( T( 129)-T( 122) ) 

UCAbt «ULAbt.*.C0( 2 56 ) • < J <1 26 > -T ( 1 23) > 

OCA bE ■ ut.Ajt.4CU ( 25b)* < T( 127) -Ti 1 2* > ) 

OC A bE ■ UC A 3 1 ♦ C 0 ! 2 36 ) * < T ( 1 2 7 ) - T ( 1 26 > > 

UCAbt-uLAbE 4.001 241 )*< T( ij2)-TI 131 ) ) 

UC AbE * utAbL + C J ( 3 1 1 > * < T < 1 59 ) - T t 1 62 > > 

UCASb>ULA^tiCul 3 1 2 > * 1 T ( 160 )-T » 163) ) 

ULAbE*ULA;»L4CD l 286) *( T < 1 5* ) -T < 15b) ) 
uCAbE*ULAjt4Co< 307 >» < T < 154 )-T ( 15b) ) 

UCAbt = OLA3t-4Co < 3 06 > * 1 I < 1 5J ) -T ( 157) ) 

OCAjE»UCA^t. + Cul 2 80)*( 1 ( Ib3)-T< 152) ) 

UC AbE* ul a:>l 4 C 0 ( 2 75)*<T(1AH)-H147)) 

OLA « UCAbt 4C0< 295 J «< T < 142 )-T( 147) ) 

UCOOL *0 • 

DO 330 I * l • 5 

UCoOL = UCOOL 4 ( CD I 2H9* I ) +CD« 34 1 ♦ 1 ) > * ( T( 14 1*1 ) -TCOOL ) 

3 JO CUNT I NUt 

UCOOL * UCOOL ♦Co< 347)*(T13l-T(141>) 

IF < NTWANb.tQ.O ) 00 TO 340 
OO 335 1 = 1 • NT nAub 
* 1 * 1 ♦ 1 O* < 1 1 4 I — iN T KAUS ) 
n2sJ29-NT«ano+ 1 
J= 1 ♦ 1 0* ( 1-1 ) 

UCOOL «uCOOL+ ( CD1M > 4CD< N2 > ) « < T < J ) -TCOOt. ) 

335 CONT I NUt 
340 CONTINUE 

UCAbE=uCAbt * ( .252/3600. ) 
vaCOOL*uCOOL« ( .252/3600. ) 

MOwLk* ( ^CAbE 4uCUOi- > >4.1 tit> 

«K| TL ( to. 350 ) UCAbt .UCOOL *MO*tW 

3b 0 F orJ.*iAT ( IrtO. • «.CAbL» • «t 10. Jt ' *C AL/bt C • * 2X . • v-COoL= • 1 10.3 « 
±> • KCAL/btC • * 

% 2x * • FOwtA = • * L 1 0 . 3 * • <w ' ) 

GO TO 2000 
4 UuO bT0P 
END 
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ju^wuC T I Nt C0-» 1 l f » 1 COWL • I caul * * CUOL *CO «L I 
CoMmoN/ I WAiMS/'l f »\A 

COMMON/ X INI T/ 1 CcM I .CUG7.CU /M.CUUV.CUVlO 
UlMLNUl ON T ( .’00 1 • 1COML ( 4 ) t CU I bOO l«c( SOO. t ) 

INITIAL TL'AP». KA To^Li»«Wt'Ol uN 1* 

I CON 1 ■ 0 

lAlc Cot<*o I co * k. • T I 

1 CON I ■ I 

UO 10 l »C»« 1 1 o * lo 
J»( I-6I/J0-M 
T ( | ) ■ TCo*t ( J I 
T ( | ♦a ) « TCAUL 

OCONbT «(T( I ) **4-T I I 4- I )••*)/( 1 ./COO 74 I ./CU7t34 l ./CllOV* 1 ./CU^l O ) 
T ( | ♦ | >■ ( T ( I ) • * 4-CCONU T /CO*>7 > •* .25 
T ( | 4-2) > I T t I ♦ 1 ) * • 4-OCONbT /CU7H 
T< | *3) a ( T I I .-2 > »*4-uC0NbT/CUbG> * * .25 
T ( | - 1 ) «T ( I ) 

T ( | -2 I «T ( I ) 

t i i - j i • r c i > 

T ( I -4 ) «T » I > 
fs. * j -^4i\ T H A NS 
IH(K.OT.IIl) oO TO 10 
TIk I »T( 1 ) 

10 CONI 1 NUL 

INITIAL TtMPtKA ToWLb. WtOlON 2 • 

DO 20 1*121,126 
T ( | ) * TC AUL 
20 CONTINUE 

T ( 1 01 ) - T CA bE 
T< 1 06 ) * T C« jF 
T » 1 Ob ) * T ( 1 1 H ) 

T l l 04 ) sT l l l Q ) 

T C j 33 » = TCAuF. 

T ( 1 02 ) * T t«jf 
T ( 127)*TCAUE 
T< | 2b) sTCAUF. 

T ( 129) ( T l 1 34 ) ♦ T CASE ) 

T( i 30 ) =0.5* < T 1 1 3 1) > ♦TCaSL > 

C INITIAL T EMpt i\ A T ONLU « WE G I ON 3. 

C 

UO 00 1 = 1 O 7 , 1 *4 1 

r ( i i a t cool 

OO luNT 1 NUl 

Jo 40 l = l««2,io3 
T ( i ) s r cajl 
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4 0 CUN T I Nut 

IF ( NTkANb.tO.O ! 00 TO 50 
0O 45 1*1 • NTKAiMO 
45 T I | ) ■ TCOOL 
50 CUNT I Nut 
kLTUKN 
UN(J 


bUBKOuT INC. HTl<T ( T * T 1M «0T tNNUbiU ) 
UlMtNblON T ( 1 JO ) »U1 200) 

Di* l I * l »NNDb 
u( I ) >0* 

I CON T I NUt 
► ETUkN 
fc.No 


SUHkCUTlNt CAP(C*T) 

0 I y t NS I UN T ( i JO ) »C ( 200 ) 
00 5 I *6 » 116*10 
C< I » »-l *0 
5 CON T I Not 

Uw c. 1*1 J7i 14. 

Cl I ) « — 1 • O* • 
el CuNT 1 Nut 

Du J l * 10* 120*10 
Cl I ) «-l .00 
.1 CuNT I Nut 

Ou 4 1 * 1 2 1 • 1 2o 

Cl I > =-l .0 
4 CUNT I NUt 

Cl 131 )*-l .00 
Cl 1 36) «-l .00 
Cl 162>*-1 *00 

wl ltoj) =- 1.00 
Cl i to 1 > =- 1 . 00 
Cl 15 b) =- 1 *00 
Cl |57)3—l . 00 
Cl 1 5? ) = - 1 .00 
Cl 142)=-1.0 
Cl 147)=-1 .0 
WtTUWN 
tND 
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SUItPOUT I ML l l.MKJtMl. 1 
CuVMON/ IpAns/nI WANS 
D I MENS I UN L ( 300 * 2 ) «LL ( bOO • 2 ) 

DATAILUJ* I ) »J*24b«264>/ 

1 l I 7* 116* i H* i 1 V* l ltt« I 1 V* 1 Jb* 1 30* 1 J4* 12SM 

2 133* 126. 1 42 ♦ 127* 126* lol • l 16* 1 lb* 1 1 J* I 1 1/ 

DATA(lL» J* 2> *J«24b*264>/ 

1 llt)*llv«i 20 *i Jb « l jo • 1 jo • 1 jO • I*. 1 • 1 2* • Ilci 

2 1 26 » 1 23 • 1 2 / » 124 t 1 j 1 • 1 j6« 1 jo « 1 jO * 1 jo • 1 Jo/ 

OfiTA(LLUiI) *J= 405*31 71/ 

1 1 62 • 1 b3 • i 34 * 1 bb • 1 b6 • 1 bto « 1 bV • 1 oO • 140* 1 b 1 • 1 b6 • 1 60 • 163/ 

DATAILLC J*2> » J * 30b *3171/ 

1 1 5 7 « 1 b 7 • i 3 d • 1 b v « 1 oO • 1 o 1 * 1 Or: « 1 03 « b • 7 » t 1 • V « lu/ 

CONDUCT JK HOUn-JPS f- ow KLO 1 ON 1 

Do 100 J> 1 • 1 i V 
L< J* 1 ) = 0 
L< J« 2 1 sJ+ 1 
luO C*>N T 1 Not 

DO 200 U= 1 20 *22* 
oo = o- i l s* 

L ( J« 1 1 > JJ 
L ( J *2 1 = 00 + 1 o 

200 CONTINUE 

«AQlATlON l>eT<*ECN SAMPLE AND HOT Ht_A T PIPE 
li-iAX* 1 2— NTH AN 3 
DO 110 1 = 1* I i-l AX 
* = 1 + 1 0* l l-l ) 

L l .< • 1 ) = N+N TpAPb* 10 
L ( K • 2 ) = N + 1 

ML CONDUCTION BtFwLEN SAMPLE AND HOT HEAT PIPE 
0=31 7+1 

L ( J « 1 > *K+NTRAnS * 10 
L< J*2)=N+1 
1 1 c CONT 1NUE 

DO 111 1=1 *12 

HAD 1 AT 1 ON LlClrtht.* HUT HtAT PIPE Ai'ju MuFFlE 
n = 4 + 10* ( I-l ) 

L ( K « 1 ) =N 
L(K«2)=K.+ l 

HE CONDUCTION bt f WEEN HOT HEAT PIPE AND MUFFLE 
*/ a 329+ I 
L ( J . I ) =N 
L ( J • 2 1 «■'>♦ 1 
1 1 l CONT INUE 

CONDuCToP HOJn— oPS FOP PEGlON 2 
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DO JOO 
JJiJ-l 10 
L ( J • 1 > ■ JJ 
L < J « 2 > ■ JJ ♦ 1 
JOO CJNT I Nub 

CONDUCTOR HOJn-UBS FO P WbO ION J 

00 400 U«265»<?90 

JJ« J-120 
L ( J • 1 1 * J J 
L ( j • 2 1 *JJ+ 1 
4UC CUNT I NUb 

uj *300 J*2vcoo4 
Jj»j- 1 a J 
L ( J » 1 > * JJ 
L< Jt2> = JJ*b 
bUO wUNT 1 NUb 

DO 700 J=24b«264 
L< J* 1 ) =LL l J • 1 ) 

L < J«2> =LL< J, 2 ) 

70o CUN T1 NUb 

L)U bOO J* J0b» 31 / 
b < J . 1 ) :LL ( J , 1 ) 

L < J • 2 > =LL ( J,2 ) 

800 CUNT 1 NUb 

I F ( NT WANS# bOi J )oJ TO bbO 
U U o<tb 1 = 1 «nFt<Ahs 

C KAq I A T l UN L>ET*EtN SAMPLE ANU CULO He AT PIPb 

K- 1 ♦ 10* ( 1 l ♦ I -.mT HANS ) 

L<K. 1 ) * < I - 1 ) * 10+ 1 
L ( < « 2 ) = 1 4 1 -N I KAN6 + I 

C Hb CONDUCTION BETWEEN SAMPLE AND COLD HEAT PIPE 

n* 32V-NTPANS+ 1 
L(K« 1 ) « ( I - t ) * 10+1 

LlK»2) s l4l— NT KAiO+ I 
m= 343+ 1 
i-< K • 1 > « 1 

H<*2)*14! -nT i<ANs+ 1 
84 b CUN T 1 Nut 
UbO CUNT 1 NUb 

DU 846 I a l . <5 

C Hb CONDUCTION BtlwbEN COLD HbAT P I Pt. ANU MUFI-LE 

J=34 1 +1 
L ( J • 1 ) = 1 J6 + 1 
L< J.2)=14» +| 

H‘«b CONTINUE 

c c^nuuc r 1 aiM bL 1 wlln hut amo culu hlat pipes 
L( 347* 1 ) = J 
L( 347.2 >= 1 41 

I TE l 6 » 600 ) » J * L l J* 1 ) »L< J» tl ) • J= 1 • J4 7 ) 
ecu FORMAT « 1 x* J I b * 7 A t J l 5 tbX« J 1 b» bA. 3 lb« bX « J I b«bX< Jl b) 
PETUWN 

End 
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bUDHOUTlNc. CO'* J ( CO • L • T ) 

COMMON/ TkAN£/ *T KANb 

CUMMON/Xll ,1 T/ ICO.^1 * CDb 7 » CD 7 b • CbBO « CUV 1 O 
JlMENbloN Cl) l 600 ) *L 1 500 • 2 ) • T 1200) 

L) I MENS I VJrx H( 10 > « XK< bOO • lb) »Kkl 10 ) 

UAT A LI • *L •* • L ‘J • C b * L 7 «LO*LO«L**LM.e-> /• AOOt # 2 * * 2 * . 32 U * » 4 00 • • JtO« 

% • '20 • • 320 * • 4 » • 2 ♦ • ‘jOO/ 

UATAl flu » 1 ) . | 3 1 . 10>/. JIG* .3 JO. .‘jbO ».OJO*1.00/*l.l fO% 1 • LbO * 2 • 7bO * 
)»J.2O0*4.'-)u0/ 

DATA Th(M1»ToIF 1*I*AL1'I MUF 1 .DEL**/ • I lb*. 010* .280 » .0 It) • 1*0/ 

c 

C AbSUMPT 1 ONb JbED IN CALCULATING CONDUCTOR VALUEb 

c 

c HtGlON 1 

C <1) .010 INCH CLEARANCE Bt. T *EEN bAMPLE AND HEAT PIPE 

C <2) .340 INCH CLtAwANCE bETwLLN HEAT PIPE AND MUFFE 1 - 

C 13) .020 INCH CLEARANCE BETWEEN MUFFLE AND HEATlNo lLEMENT 

C 1 4 ) .480 INCH CllAPANCE PLTrfELN HEATING ELEM + NT «-D FIWST bTlFFNE*. 

C lb) . 010 STIFF i*LH THlCKNEbb 

C 

C 

C • NbM iiNJJ Ai<E I HE NUMBEk OF L AYtki OF lUNGbfEN »MOLY AND bT bTEE 

C 

blG/!A= 1.71 4E- 09 
Nbw 3 30 
NbM=30 
NbS*30 
C 

C I HE UN 1 • J UN i* ♦ DL.ua I H 1 P. * 1 b I F * * M Ai_ • 1 I»>UF Ape F I • • He. V AlUE-> Ape 

C pLAD IN IN I NCrltb AND CONVERTED TO FI. 

C DElX lb AXIAL DISTANCE STEP b I EE 

C T MuF lb THE MUFFLc. THlCKNEbb 

C ThTW lb THE .3L A T EH THICKNESS 

C TbTF lb THE STIFFEN THlCKNEbb 

C 

c 


c 

table no. 

MATEP 1 AL 

PKOPEP 

c 




c 

i 

bT. STEEL 

lP 

c 

2 

bT. STEEL 

K 

c 

3 

AlUi'i 1 NA 

CP 

c 

4 

ALUM 1 N A 

K 

c 

■J 

TUNUbTEN 

CP 

c 

o 

T UNGbTEN 

K 

c 

7 

••iOLY 

CP 

c 

8 

MOEY 

K 

c 

•y 

T M-U/W 

UP 

c 

10 

T A — kj/ w 

K 

c 

i 1 

SAMPLE 

CP 

c 

IP. 

sample 

P 

c 

13 

HEAT PIPE 

CP 

c 

14 

HEAT PIPE 

K 

c 

16 

COLD H.P. 

CP 

c 

16 

COLD H.P. 

K 
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X<vib = 0 » 06 J 
lCmP* .400 
DO 10 I * l ' l 0 


10 K ( I >*«w< 1 )/| 2 « 

DbLX*UELXX/| 2 v 
THTW*THTKl /\ 2 • 

TbTF*TbTFl/i2. 

T«AL*T wAl 1 / 1 2 • 

T i»hJK * T >*t OF 1 / i 2 • 

Pi * J, 14ljV27 
A2*2.*P1 *P< ? ) *DtLX 
AJ*2» *P i *P t 3 ) *0tLX 
, A4 = 2»«Pl»P(4)’*0toX 
Ab*2» *Pi *k l 5 ) *0tLX 
A6=2.*PI *W<6) «L)LLX 
A7s 2 • *P l *P ( 7 ) *L)LLX 
Att = 2. *P1 «P ( ) * DLlX 
AV = 2.*P1 *1 (Q ) *UtLX 
A 1 0 *2 • *P l *P < 10) *ytLX 
Hl2b = b.*UELA 
*126 = J«b*DtLX 
*131 = * ( 7 ) 

AN2Jb = PlMK12y«-P126)*DLLX 

AM2b* = PI * ( K120+P1 31 >*ObLX 
F2l »K( 1 )/*<2 ) 

Fa^skUl /* ( b ) 

F6«5sW(b)/P(fi) 

FI l |*1 .0 
FI 1 b= 1 • C 
FI 10*1 .0 
F 127* 1 . w 
F 1 16*1 *0 
F 1 l J= 1 • 0 
FI IV* 1 .0 
F 1 20= 1 • o 
C 

c CALL iNTfcWPOLAl ION *UUTlNb 1 Im P*tVlOUb T LMPL*A T U*db TO Lb T tkM 1 NL 

C Nbw CONDUCTIVITY 

C 

OO 700 J = 2 • 1 6 * 2 
DO 700 1=1 .317 
N1 =L ( I . 1 ) 

\2 = L< 1 .2) 

I mvG= ( T INI ) * I * N2 ) ) /2» 0 
C Aj_L 1NTPI TaV0» J*F \K> 

Xiv ( l • .1 ) =XivK 
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KiO CONTINOL 

copi « < i •/ i i • /Li* i */c i -i • n *:»igma*A 2 *F 2 1 

CD 12 » AM 1 , 14 ) *PI MW< J)+W< 2 > )/»H< J)-w»2 > >«DLCX 

LL«|J « AM 1,14) • K| *< W< 4 ) +W ( J ) ) / < Wl 4 ) -W < .1 > > •DtCX 

CD«S4M t .✓» I , / L •> + 1 • /L4- 1 • ) > «bloMA«A*j<»Fb4 

CDfc'j « < 1 ./( 1 */L6+l ./Lb-1 • > ) «bl OMA«A6*Kt»a 

CDfi 7 * < 1 •/ 1 1 */L6+ 1 */L7-l • ) > *6lGMA«A6 

CD7e * <C*/I 2*-t* ) )/» NSKf-M 2. *fcW ) /C2*-tW )*< 1 */fc7+ | */bW- I • > >* 

» 3 1 Gma * A 7 

CDQ9 ■ i£'*/ < 2 *-trt) )/l NSM+4 2**tM)/<2*-EM) M l */t8+l */EM-l * ) > * 

» jMA*Ao 

CD910 « 1 t3/ I 2 • ”LS» ) >/ <N6i + < 2* *tS)/ < <!.-Li ) * < 1 */t9+l ./£a- 1 • ) ) * 

A ;>IGmA*A9 

I F ( I C of 4 1 (tWtO)KbT OWN 
A237«PI * < W ( | 0 > * *2“W ( 7 > **2 ) 


LL>2 j7sCu7tt/A /»A^J7<P 12U 

CD26 1 » IU/I 1 ./£6 */ES> > *SIG*A»PI 

CD290 * < 1 •/ l 1 • /E5 + I ,/ECHP ) ) *3l GMA*A5*KbA 

CD296 ’ < I • / < 1 • /ES+ l */E7) ) *S IG* V 'A»A7*R( 5 )/W< 7 ) 

Cl) jO J « CU7h*A 8/A /*w( 7)/w( t>) 

COJ06 « CDG r ) + A9 / Ao+k < ti ) /w ( V ) 

CD J 1 1 = ll)V 1 j»A10/Asy*W(S»)/w( JO) 

Ow 1 00 1 * l * J i i * 1 0 
>v- I ♦ I0» <1-1 4iM T ) 

1)- < a. *Ct . 1 1 1 ) 

*CD< I ) *CU2 1 *< F 1 k ) **2 + T < 3 ) * # 2 > * < 1 ( k ) + T ( 3 ) ) 

CD < I + 1 ) *C0 J2 * A k 1 1 + 1 • 1 4 )/XM 1 • 14 ) 

CO< l +2)=C043*Akl 1 +2 • 14) /XK < 1*14) 

CD ( 1 ♦3 ) =C054 # * T 1 I +3 ) + +2 + T < 1+4)**2)*<T< 1 +3 ) + T< 1+4) ) 

CD < 1 WD sCL)6 r i * l T ( I +4 ) * * 2 +T < 1 +“ ) **2 ) M T < 1 +4 ) +T < I+*j) ) 

CD < I + b ) »CJ67* ( T ( 1 +5 ) * *2 + T < 1 +6) **2 > < T< I +b > + T < 1+6) ) 
CU< I +6 ) * CD 70 * < T < I +6 ) * *2+7 < 1+7)**2>*<T< | 46 ) + T ( H7) ) 

CD < 1 + 7 ) sCJdo* < T( I + 7 ) v *2 + T < I+o)**2)*<T< l+7)+T< 1+tJ) > 

CD ( 1 4fa ) = CD9 i J * { I 1 1 -*-0 ) **2 + T < l+y)**2)*<T< 1+6 )t 1 < 1+9) ) 
100 CONTI Not 

CD ( 1 20) =AK < 1 20 » 12 ) *PI *W< 1 ) **2'DECX 
oO( 121 )=0*0 

CD ( 122)= AM 1 22 » 14 > *P I * < P< 4 ) **2-W < 3) **2 > /DtCX 
CD < 1 23 ) =0 • 0 

CL)< 124)=AA<i24«iO)*P|*<(W<b) + TNiUF ) *• 2-W < b ) * *2 )/DELX 
CD ( 1 25 ) =XK ( l2tl»4)«P|*<<W<6> + THTk) *«2-k< 6 ) **2 >/CELX 
CD ( 1 26 ) = XI < j 26 » 10 ) *PI * < < W< 7 ) + TETF ) *»2-W< 7 ) **2 >/DELX 
CD < 127)= AM l 2 7 * 1 0 ) *P 1 M t tt < Q ) + T b T F ) **2-K< G) * + 2 )/DtLX 
CO I 1 20 ) =XK < i 2 8 • I J ) *PI * ( <k< 9) +TS.TF ) **2“k< 9 > **2 ) /DtCX 
CD ( 1 29 ) = AM i 29 ♦ 2 ) *P I * < < W < 1 0 ) + 1 d»AL ) **2-w < l 0 ) **2 ) /OELX 
DO 200 1 « 1 20 t C 1 J * 1 0 

U)( 1 + 10) SCO < 120) *XK< l+10» 12)/AK( 120* 12) 

CD ( 1+11)80*0 

CD ( l+12)a CD <lc2) +XK. <l + 12«l4)/AK(lt2*l4) 
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CO ( I ♦ 1 j > * o . o 
CO< I ♦ 14 ) *CO( 1 24 ) »XK I 
CO ( l 4 lb) »CD( 1 2 b ) *XK ( 
CD ( |416)*CO( 126 ) *XK( 
CO( 1 4 1 7 ) «CU( 1^7 ) »XK ( 
CO< I 4 lU)>CO( I2tt ) *XK ( 
CD< I 4 J9 I «CO( 1^ ) *XK ( 
<!00 CONTI NUL 


|4J4* 10)/XK( 124* 10) 
1 4 | b*4 ) /XK ( 125*4 ) 

I 4 16* )0)/XN( 126* 10) 
1417* I0)/xk( 127* 10) 
1410* I0)/Xk( 120* 10) 
1 4 1 V • 2 >/XK( 120.2) 


CoNUeCTuW VALOLO Ko H ktGlON 2 


CJ( 200) 
C0( 231 > * 
C0( 232 ) * 
CDl 233) = 
CO(2J4 ) = 
CD( 235) 
CDl 2 36 > = 
CD< 237) 
Co< 230 ) 
CDl 2.'9 ) 
CD< 240 ) 
Cl)( 24 1 > 
CU( 242 ) 


« x« (230*2) *P l« ( (k( 10 ) 4T*AL )*«2-W( 1 0 ) **2 > /DtCX 
CD( 2 30 ) *Xk( 2 ji .2) /XK( 2 30. 2) 

CO ( 2 30) •Xk(232»2)/’XK(230.2) 

CO ( 2 30 ) *XN < 2 33* 2 ) /XK ( 230 « 2 ) 

C0( 230) *XK (234 »2>/XK< 230.2) 

» X<(233»2)»AM23b/(W125-K126) 

CU2 3 7 * « T ( 1 26 ) * * 2+ T ( 127)**2)*(T( 126) 4T ( 127) ) 

» CD237* ( T( l2ti)**24T( 127)**2)*( T( 128>4T( ,27) ) 

* C02 3 7* l T( 129 ) **24T ( 12rt ) **2 ) * ( T ( 120)4 T ( 1 28 ) > 
a Co2 3 7* ( T ( 130) *«24T ( 12V)**2)*(T( 130)4 T( 129) ) 

= 0 • J 

a CO ( 126 )/2.*AK( 24 1 . 10)/ XK ( 126*10) 

= CDl 126)/2.»XK(242* 10)/XK( 126*10) 


Cl) ( 243) 
CD( 244 ) 
CO ( 245) 
CO ( 246) 
C0( 247 ) 
CO( 248 ) 
C_M 249 ) 
CO( 250 ) 
CO( *:51 ) 
CO ( 252 ) 
CD( 253 ) 
CO( 254 t 
l.O( 255) 
CO< 256) 
«-0( 257 > 

CD( 258 ) 
CD( 259) 
CO ( 260 ) 
CD( 261 ) 

1 


* CDl 126 )/2.»x*(243» 10 )/Xk( 126*10) 

3 CDl 126 ) /2 • ♦ XK (244* 1 0 ) /Xk ( 126*10) 
a x<( 24 5 » 1 0 ) *P I » ( W ( 8 > 4k ( 7 > )*T5TF/<W<8)-k( 7) ) 
a X< (246*10) *P 1 * ( k ( 9 ) 4W ( G ) ) *TSTF/ ( k ( 9 )-k( 8) ) 
a Xk ( 247 • 10)4 PI*(k( l 0 > 4W ( 9 ) ) * r5TF/(R( l0)-k( 9) ) 

= CDl 2 30 ) »XK.( £>40*2 ) /Xk ( 230*2 ) 

= CU23 7/2* 1‘118) **24 T ( i 30 ) **2 ) « l T ( 1 18)4 T ( 130 ) ) »F 1 18/F 1 ?o 

- CD2 37/2 * # ‘T ( j l 9) ** 24 1 l 1 30 > *»2 > * ‘ T < l 1 9 ) 4 r l 130 ) )*F1 l 9/F 1 28 
a XK ( 251 » 10)*PI# (K| 264). I 7) ) *T 3TF / (kl26-k( 7) ) 

= <K < 252*10)i*t->l*(«( lC)4kl26)*T5TF/(K( 10) -k 126) 
a CDl 2b 1 ) *Xk l 253 » 1 0 ) /xk ( 251 * 1U> 
a CDl 2 b 2 ) *Xk ( 254 * 1 0 ) /Xk ( 252 • 10) 

= CDl 251 ) »Xk(255* 10)/XK(251 » 10) 

= CD l 2b2 ) * XK l 256* 1 0 ) /xk ( 282 » 1 0 ) 

= CDl 25 1 ) *XK( 257* 10) /XK( 251 » 10) 

= CDl 221 )*XK(258» 10 )/Xk( 252» 10) 
a XK ( 259»2)»AM259/lR126-kl31 ) 
a XK l 260 *2 )*p 1 *k 1 3 1 « «2*DECX/W 131 
a CD 2 6 1*1 IWI6)4 ThTW)» * 2-k ( 6 ) **2 ) *F 1 16* 

' ( T ( 1 16 ) **24T ( 1 36)**2)« ( T ( 1 16)4T ( l J6) ) 


CO ( 2 62 ) = C026 1 * 1 1 k ( 5 ) .-TMvjF ) **2-k ( b > **2 > *f- l l 5* 

1 l T ( 1 i p ) * *2+ T ( 1 36) *«2) *1 T ( 1 lb)4T ( 136) ) 


Cu(26J) = CD 2 6 1* Ik( 4) **2 _h 'l?)* # 2) # Fl 13* 

1 (T ( 1 13)**24l ( 136)**2)»( T( l 13)*T( 136) ) 

CO ( 264 ) = CD26 1 *l*< 1 ) **2*F 1 1 l * 

1 ( T ( 1 1 1 ) »*24T ( 1 36)**2) * l Tl I 1 1 )*T( 136) ) 

IF ( NT PANS* GT • 0 ) CO <264 )~0» 
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c 

C CUNUUCTOWj FOK ULOION 3 

C 

CUC t'Ob) *X*. ( ;»(,a » ib ) «H| • ( W ( 4 ) ♦ »< « 3 ) )/UL«-X 

cuuob ) =cu < ?o*j ) * Xn ( 2bb» 1 1> >/xk e 2bb» 16 > 

CU1 t'67)»Cul?b3l *xx»i , c>7« lb)/XKe 26 b» 16) 
cue cCti ) =cue 26 b ) *Xk < 260. 16 > /xk e 205 * 10 > 
cue cbs» ) =o*o 

cue 270 > *cue i 24 ) *xk< 27u» io>/xk e 124 • 10 ) 
cue 271 > «coe 1 2 4 ) »xk e 271 • 1 o>/xk e 124 . 10 > 
cue 272) = CU< 1 24 ) *XK( 272 » 1 0>/XK< 124.10) 
cue 2 73 > =CI e 1 24 ) *XK(273. 10 ) /XKe 124.10) 
cue 274 >*0.0 

CDe 27b) *cue 1 2 b ) * XK e 27b. lOl/XM 126*10) 

CUe27b)*CU( 1 2 6 ) * Xk e 27b . 10 )/Xk< 12b. 10) 
cue 277) *cue 1 20 ) *XK< 277. 10)/Xxe 12b. 10) 
cue c 76)=cue i 2o ) *XKe27t). i o > /xk e 1 2b . 10 ) 
cue 27v) =0 • 0 

v.ue 2M0) =cu( i 2 / ) *xxe 2b0. io >/xk( 127. io> 
cue 201 ) =cue i 2 7 ) *xk< 201 . 1 o>/xk e 127. 10 ) 
cue 2H2 > =cu( 1 2 7 ) *xx e 202* iu>/xxe 127*10) 
cue 203 ) =cue 1 2 7 ) * Xk ( 203. 1 o> /xxe 1 27 . 10 ) 
cue 204 > =o#o 

cue 20b) * 1 . 0O*AK.< 20 b . 2 )*PI *UELX 
cue 206 ) =cue 1 28 ) *xk e ?eo. io>/XKe 120.10) 
cue 2 B 7 > =co< 1 2 b ) *xxe 207. 1 o> /xk e 120*10) 
cue 200> =0.0 

cue 20 V) =cu e 1 2 9 ) *xxe 200.2 )/xk t 12V.2) 

COe 2 VO) = CD290* ( T e 137) **2+T « 142) *«2)* e T e 137) ♦H 142 ) ) 
CUe 2 Vl) = CD290* l T e 1 30 ) ** 2 + 7 e 143 ) ** 2 ) * e T e 1 30 > + 7 e 143 ) ) 
CU12V2) = CO290* 1 T e 13V)**2 + 1 < 144>**2)*< 1 ‘ 13V > + 7 e 144 ) ) 

CU e 2 V 3 ) * CU 2 VO* 1 1 e 140 ) ** 2 + J e 1 45 ) ** 2 ) * < T < 1 40 ) + T e 1 45 ) ) 
cue294) = CD290* ( T e 1 4 1 ) ** 2 +T e l 4 b ) ** 2 ) * < T « 1 4 l ) + T e 1 46 ) ) 
cue p o** ) «x< e ^95 » 2 ) *pj # e we 7 ) +«eb ) ) *delx* . 50 / e w e 7 ) +w e 5 > > 
cue296) = CD296* (Te 143)«»2+Te 148 )** 2 )*< T< 143)+Te 148) ) 
CU < 2 V 7 ) = CD296* ( Te 144)**2 + T< 149)**2)*( T e 144 )+T( 149) ) 
cue290) = CD296* ( T ( 1 4b) **2+T • 1 50 ) * * 2 ) * e T e 1 45 ) + T e 1 so ) ) 
CD < 299 ) = CD296* (Te 146)**2 + T e lbl )** 2 )*e T ( 146> + Te 151 ) ) 
cue 300 )=CUei>9b)*2»C*XKe 300.2 )/XK( 2 V 5 « 2 ) 
cuejon = CD30 1 * 1 T e 140) **2 +t e is3 ) *# 2 >* e t e i4b) + T e ib3) ) 
CU e 302 ) = CD3 j 1* i T ( 149)*«2 + T e 154 ) *#2 ) * e 7 e 149 ) + T e 154 ) ) 
CU1303) = CD30 1 * t 7 e 1 50 ) ** 2 + T e 1 55 ) *# 2 ) * < T e 1 5 C ) +T e l 55 ) j 
cuej 04 ) = CD301* ( T< 151 )**2+T ( 156>*»2)*( T ( 151 )+T( 156 ) ) 
cue J05 ) =cu ( 28b ) * xk e 305. 2 >/xk e 2 b 5 * 2 ) 

cue:i06) = CD3 J6* < 7 e 153) **2 + 7 e 1 57)**2 )* e 7 e 153 ) -*-7 < ) 57 ) ) 
cue307) = CD706 * 1 T e 1 54 ) ** 2+7 e l 50 ) ** 2 ) * ( 7 t l b4 ) +7 e 1 50 ) ) 
CUe30H) = Cd?06* 1 7 ( 1 55) **2 + 7 e l 59 > **2 ) * e 7 e 1 55 ) +7 ; 1 59 ) ) 
cue309) = CD3 06* 1 7 e 150 ) *+2+7 e l 60 ) **2 ) * e T e 1 5b ) +7 ( 1 60 ) ) 
cue3ioi=xK( 7iJ»2)*ui*(R < 1 0 > +w e 9 > > *uelx/ ewe 1 0 ) +we 9 > ) 
cueuil) : CD3 1 i * ( 7 e i59)**2+7 ( 162 >**2 )* e T ( 159 )+7 e 162 ) ) 
cue312) = cr>3 1 1 * ( 7 e 1 60 ) **2 + 7 ( 163 ) **2 )* < 7 e 160 ) + 7 e 1 63 ) ) 
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CUt 31 31 = CO t J 24 1 »XM 31 3. 10) /XK ( 124.101 
Cot J 14 ) »Cu t 1 2 b 1 »XN ( J1 4 % I O ) /X* t 126* 1U I 
cut .11 b> «Cu t 1 2 7 1 * AK l 31 b» I0)/XKI 127. 10 I 
CUt JI6)*CDt 1 2b I »AM Jib. lUI/XKl 1 £b« 10) 

Cut Jl 71 =CD t i 2S> 1 • A*v‘ 31 7.2 1/XM12V.2 ) 

CDt 2611=0.0 

cut 262) >0.0 
cut 26J) >0.0 

CUt 264 1=0.0 

IF | NTWANb.LO .0 1 GO TO 2b0 
L)vJ 24b 1 = 1 . N T ixA'Xb 
*■ 1 ♦ 10* I 1 1 ♦ I -nTNANb 1 
j« 1 ♦ 10* t l - I I 

cut M «CU2l * t T ( J 1 **2 + T t 137 I ** 21 * t T t J| +T t 1 J 71 J 
24 b uONTlNOL 
2b0 CONTINoL 

IMF F uluj. 1 NG CUiMUOCT UHb A«L Mjk 1 ML HELIUM GAb CON T A I NLU In 
TMF MUFFi-E 
Al«2.*Pl»Ml J )*DtLX 

AVGl *t Al*A2>/2. 

XU|S1«W<2)-R( 1 1 
AVG2* t A4-.AL ) / 2 . 

XU! ii?*Wlbl-Rl 4) 

CON ML I =XNMt * A Vo 1 / Xu 1 b 1 
CONML2* XilML # A Vo 2 / XU 1 b2 
UO 2001 1= 32V 

CUt 1 1 *COWHt 1 

2001 oONTlNUL 

UO 2002 1«J30»346 
CUt 1 1 »CCNHL? 

2002 CONTI NOt 

C CONUOCTION BElwttN MOT AND COLD HLAT MIPLb 

NCAYEW= JO 
LC AYLM = Lm1 

CUt j4 7)=PI*tKt4)**2-N‘2)**2>*i»l GMA* t LCAYLM/ t 2 .-LL AYLP 1 1 /NLAYLP* 
i* < T t 3 ) **2 + T | 141 1 * *2 1 * ( T l J 1 + T t 141 1 1 
mLTUWN 
LND 
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-»OfiHQOT 1 NE To I D l C • v»D « NC T OL • M (NNJj * NMH • Mh |N«T l i*i * l PT * T bCON • UT ) 

U I wLNj I ON C< £00) ♦CUCbOO) • M ( £00 • lb) 

bWAL a I • QE30 

DO 100 I ■ 1 tl«ML)S 

CAP « C< 1 ) 

IK(CAP) 1 00 • 1 J 0 • COO 

400 continue 

NCONO « M ( I « 1 ) 

Cbu'*' >0*0 
DO £00 o ■ 1 •rtCO.NU 

1C » J*£+l 
Nc a M ( I « l C ) 

CbUM * CbU^ ♦ CD (NO 

£0O CONTINUE 

bTEbT a CAP/CoD>» 

IF< STEbf .GL.O'tAL ) GO To 300 
bMAC = STEbT 
NCT*I_ a I 
300 CONTINUE 
IOu CONTINUE 

OT a SMAc/TsCoN 

TOP a f- L'JA T ( | N 

TT 1 a TIM+OT 

IDT a 0 

IF ( TPW.GT . TT i ) oO TO bOO 
UT a TPK-T I Ml 

IOT a 1 

buO CUNTINUE 
WETUPN 
ENo 


A - 1 7 


LOCKHEED HUNTSVILLE RESEARCH A ENGINEERING CENTER 



LMSC-HREC TR D697553 


bURWOUT I Nt |NTP( X • N • Y ) 

Cumnon/uT A/ cC l 40 »bl ) 

u m CC ( N • 1 ) 

Jl a J** 

Xk-ST ■ CC(N,2) 

XHST • CC(N'Ji) 

IF | X.EO. XL ST > GO TO 700 
1F< X.lT .Xt_6T 1 GO TO 100 
IF ( X.GT.XMbT ) GO TO 200 
JT » 4 

buO CUNT I NUt 

I F ( CC ( N » JT ) — X ) 300*600.400 
JOO CUNT I NUL 
JT ■ JT +2 
GO TO 600 
400 CUNTlNUt 
N 1 « JT-2 
N2 « N 1 ♦ 1 
NJ ■ JT 
N4 = N J+ 1 
GO TO 900 
6o0 GONTINUt 
NT a JT ♦ 1 

Y = CC(N.NT) 

GO TO 1000 

700 CONTI NUt 

Y a CC ( N« 3 ) 

GO TO 1000 

ioo continue 
ni a «> 

N2 * J 
N J - 4 

N4 * b 
GU TO 900 
2u0 CUNTlNUt 

Nl = 2* J-2 
N2 * Nl + 1 
NJ = 2* J 
N4 a N3>1 
900 CONT 1 NUt 

uo a (CC(N«N4)-CC‘N.N2 ) )/(CC(N.nJ)-CC(N«N 1 ) ) uu 

Y = CC1N.N2 ) +JL* (X-CC(N»N1 ) ) 

1000 CUNTlNUt 

HETUWN 

END 
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bOHWOOTINE DATA 
COMMON/oTA/CO <* 0t51 » 

701 FOHMATI2IS) 

702 T OOMAT ( 2U 1 0.0 ) 

700 KUHWAT < 1H| « | JA * i IHI A h L t b) 
7o<* KOHMATI lHO*bX*2l7) 

VOt> FORMAT ( IHO »sx »2t |2»S> 

W«l TL (6* 70 J) 

100 CONT I Not 

KtAOt ta« 701 )N» J 
| T6 ( t>* 704 )N» J 
IFcN.EO.Ol GO TO 000 
CC(N*li « J 
DO 200 1 ■ 1 * J 
HtA0<S»702> A » 0 
*W!TE<6»7O5)A»0 
is 1 « 1*2 

ls2 « IS 1 ♦ l 
CC ( N « K 1 ) > A 
CC(N.K2> * B 
200 CONTI NOt 
GO TO 100 
•) CONTI NOL 
WtTUWN 
ENo 
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too* uo | JO 


0*0 too* 


20 Ou • 200u* 200 0* 2000* 


b 


1 

V 


0*0 


• Jbb 

200 • 


• oob 

sue. 


• 100 

7o0. 


• 123 

1000* 


• 133 

1 bwO • 


• 1 A2 

2 CoO* 


• 1 b2 

2bwC • 


• IbA 

3000* 


• 177 

2 

V 


0*0 


A.b 

2a0* 


0.7 

bOO* 


0*3 

7 j0 • 


v*o 

1000* 


11*! 

IbUO* 


13*2 

2000* 


1 b*2 

2b00. 


1 7* 1 

30u0. 


l->*0 

3 

10 


0*0 


* 0 1 6 

200 • 


* 022 

ACC • 


*025 

ooo* 


• 027 

bOO* 


*028 

1 COO* 


• 02b 

1 bOO • 


• 0 30 

2000* 


• 031 

2b00 • 


• 032 

3000* 


• 032 

1 J 

2 


1 00* 


* 100 

AOU0* 


• 100 

1 A 

4m 


100* 


1 *0 

A CoO • 


1 .0 

1 b 

(1 


100* 


• 100 

A 000 • 


• 100 

lb 

d 


ICO. 


1 *0 

A 000 • 


1 • 0 


CP - ST A 1 NLt 3b bTfctL 
< - STAlNLtbb STfctL 
CP - ALUMINA 


CP-MtAT PiPt 
iC-HKAT PI Pd 
CP-COLL) HtAT PI PL 
k--CG^U HtAT PIPt 


'j<*0 • 


•jAO. 
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4 

10 




0*0 


22. b 

A ~ 

Alum 1 NA 

tuOt 


1 /.6 



400. 


1 < . 6 



600* 


V. J 



oGO » 


/ • 2 



1000* 


b. 7 



1‘juO. 


3. V 



2 OGO. 


J.2 



2bO0. 


3.2 



3000. 


3.1 



b 

2 




i oo. 


.036 

CM 

- TuNObTLN 

3bo0 * 


.036 



6 

10 




luo. 


22o.O 

K - 

T UNObTLN 

2 oO» 


1 30 • O 



300. 


10b. 



bus. 


93. 



1 OuO • 


7*. 



1 bOO • 


73. 



2000* 


6b. 



2bOO . 


6b. 



3000. 


62. 



JbOG » 


b9.b 



7 

12 




0- 


• 046 

OH 

- .HOLY 

C'JO • 


. 0b7 



bOO . 


.062 



7bO • 


.064 



1 G'O. 


• 0 65 



1 2bC • 


. 0 66 



lbOO. 


. 067 



2w '0* 


.070 



2500 • 


.0 76 



3000. 


. 063 



3bu0 • 


. 091 



4 OuO ■ 


• 100 



a 

1 2 




0.0 


64 . 

N - 

■ MULY 

2b 0 • 


62 • b 



bOO. 


7-*. 



7b0 • 


7o. b 



J OilO. 


V 4 • b 



1 2 b 0 * 


72. b 



1 boO* 


ov. b 



20UC. 


63 . b 



2bOC. 


b 7 , 



3000 * 


bl .0 



JbOO * 


4b. b 



4000. 


42. 
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V 

2 




0*0 


• 1 00 

OR* - 

T A- 10* 



• 1O0 



lo 

b 




0*0 


Jj»b 

K - 

T A- 10* 

(.'400* 


JJ.b 



J.’OO* 


2V«b 



•4 000 . 


2o« 



•*U00 • 


22 «b 



1 1 

J 

4- 




100. 


• O Jfr 

CM - 

TUNoblLNi SAMRLL ) 

JbOO. 


• 0 Jft 



1 2 

1 o 




i oo* 


220.0 

s - 

1 uNGb T t_N 

joo. 


1 JO. 0 



000. 


l 0 ->. 



‘>uO . 


O J. 



1 ooo • 


/V. 



1 boo. 


/ 



2000. 


t>ti • 



2bO0 • 


o.>« 



JOOO. 


b<- • 



JbOO . 


t>‘V • 
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